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Abstract 

We discuss the observable effects of enhanced black-hole mass loss in a black hole- 
neutron star (BH-NS) binary, due to the presence of a warped extra spatial dimension of 
curvature radius L in the braneworld scenario. For some masses and orbital parameters in 
the expected ranges the binary components would outspiral, the opposite of the behavior 
due to energy loss from gravitational radiation alone. If the NS is a pulsar, observations 
of the rate of change of the orbital period with a precision obtained for the Binary Pulsar 
B1913+16 could easily detect the effect of mass loss. For M BH = 7M & , M NS = 1.4M , 
eccentricity e = 0.1, and L = Wum, the critical orbital period dividing systems which 
inspiral from systems which outspiral is P~6.5 hours, which is within the range of expected 
orbital periods; this value drops to P~4.2 hours for Mbh = 5M . Observations of a BH 
pulsar system could set considerably better limits on L in these braneworld models than 
could be determined by torsion-balance gravity experiments in the foreseeable future. 
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1 Introduction 



Extra spatial dimensions, beyond the three encountered in everyday experience, have long been 
discussed in theoretical physics [HI HI [29] . Earth-based tests of these ideas are difficult [TJ. We 
discuss an astrophysical test for the consequences of a particular class of extra spatial dimension 
models. The observations and analysis would be similar to work done on the Binary Pulsar 
PSR B1913+16 (two neutron stars, one observed as a pulsar), which provided a high-precision 
test of the effect of energy loss by gravitational radiation [3TJ [321 [33] . In the test presented 
here the binary pair would consist of a black hole (BH) and a neutron star (NS) observed 
as a pulsar. The intended audience for this paper comprises both astronomers familiar with 
observational astrophysics but unfamiliar with extra dimensions, and theoretical physicists in 
the opposite predicament. We hope that this paper and the similar works which have proceeded 
it [61 [27J [121 H31 [22] will help these two groups to bridge the gap between them, and demonstrate 
that there are potential astrophysical avenues for progress in this area. 

One reason for interest in extra dimensions is string theory, the most notable quantum 
gravity theory, which requires 6 or 7 extra spatial dimensions [25] . Such extra dimensions could 
evade detection if they are extremely small — "compactified" or rolled-up; a coordinate r along 
such a dimension is periodic on a length scale L, i.e., r = r + L. The most natural scale for 
these dimensions is the Planck length L ~ a/ HG / c 3 ~ 10~ 33 cm. Probing that length scale in 
accelerators would require energies at the Planck energy ~ ^[hxfjG ~ 10 19 GeV, 10 16 times 
larger than the electroweak TeV scale of current experiments [5] . 

A second motivation for considering extra spatial dimensions is the "hierarchy problem" - 
the relative weakness of gravity compared to the other three fundamental forces. If gravitons 
alone propagate in extra spatial dimensions, then the gravitational field of a particle drops faster 
than 1 /r 2 on length scales smaller than the size of an extra dimension, and is therefore weaker 
on larger scales where the inverse- square- law behavior becomes manifest. In addition, L could 
be ~ 1 mm and be undetected in particle accelerator experiments (which don't probe gravity), 
and torsion-balance gravity experiments [3j [2J . Gravitational torsion balance experiments have 
set a limit of L < 44//m (95% confidence) [15]. [1] state that modest improvements in these 
torsion balance experiments will occur, but their own projections show sensitivity far beyond 
the current limits may prove difficult in such experiments. 

An alternative to compactification is possible. Extra dimensions of infinite size are allowable, 
if their geometry is "warped" so gravity cannot propagate further than the length scale set by 
the torsion-balance experiments. For one class of models all forces and particles other than 
gravitons exist only on a "brane," which is the boundary of a infinite, warped "bulk" (e.g., 
the Randall- Sundrum 2 model, RS2) [28J. In this "braneworld" scenario, our 3+1 dimensional 
world is the brane, populated with the standard-model particles and forces. An extra spatial 
dimension is "orthogonal" to the brane, forming the bulk. In such warped models the bulk 
must be an anti de Sitter (AdS) space. Gravitons propagate on the brane and in the bulk, 
thus, in principle, resolving the hierarchy problem. The "warping" amounts to a redshifting 
of gravity as it propagates into the bulk, so, in effect, gravitational experiments only probe a 
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short distance into the bulk. 

Much work has been done to understand the nature of black holes in the braneworld scenario 
[9]. Because the bulk space is anti de Sitter space, application of the AdS/CFT conjecture [2"U] 
becomes possible, particularly in the context of RS2. The result is that the full classical 5D 
braneworld black hole solution is equivalent to a 4D quantum corrected black hole. Moreover, 
this analysis yields a dramatically increased evaporation rate for large black holes due to the 
existence of conformal degrees of freedom [7j. This conclusion was challenged by [S] in which 
a static black hole solution was demonstrated to exist. This solution however was found to 
be unstable. Moreover, numerical investigations have failed to yield evidence of a stable static 
macroscopic black hole in the context of the braneworld scenario [191 El] • The lack of a dynam- 
ically stable solution may serve to indicate that the RS2 model implies a dramatic increase in 
the amount of Hawking radiation emitted by macroscopic black holes. 

Note that we can give a broader heuristic argument for the general significance of quantum 
gravity in relativistic astrophysical situations involving black holes, without tying ourselves 
directly to specific models involving extra dimensions. On general grounds [161 E], quantum 
gravity, due to its essential non-locality dictated by the underlying classical symmetry of general 
coordinate transformations, can be argued to be sensitive to other scales, such as the TeV scale 
of particle physics, where, also, by unitarity, new physical phenomena are expected to appear. 
Similarly, the increased evaporation rate can be expected based on models of black holes as 
bound states of these essentially non-local quantum gravitational degrees of freedom. Such 
models often arise in the analysis of black hole evaporation in string theory [2T], in which the 
non-local, stringy, quantum gravitational degrees of freedom can spread from the Planck scale 
to large, astrophysically significant scales. Thus a very large number of degrees of freedom can 
be transferred from very small to very large length scales, providing an enhanced evaporation 
rate that is crucial for the discussion that follows. 

2 Astronomical Limits 

Some observable astrophysical implications of these results were explored by [5J. Observations 
of X-ray binary systems (BH + companion star) can potentially constrain the size of extra di- 
mensions in the braneworld scenario. The constraints rely on the kinematic age for the system 
or observational limits on any change in the orbital period of the system. These observations 
have the potential to set lower limits than can be obtained by current torsion-balance exper- 
iments. Using the kinematic age of the X-ray binary XTE J1118+480 [27J set an upper limit 
of L < 80/im (95% confidence). [12] set an upper limit of L < 970/im (3a) from current ob- 
servational limits on the change in orbital period of XTE J1118+480, while [13] set a limit 
of L < 161/xm (3a) from observations of A0620— 00. [12] states that one additional measure- 
ment of the orbital period of XTE J1118+480 could constrain L to less than 35yum. However, 
the physical complexity of such systems and our limited understanding of their astrophysical 
behavior makes a proper interpretation of the observations difficult, and perhaps ambiguous. 
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While torsion-balance experiments may not be able to set such low limits, they provide cleaner 
experimental setups, with a more straight-forward interpretation. 

Cleaner astrophysical binary systems comprise components that can be treated as point 
masses, with no mass exchange. [22] discusses a BH BH or BH-NS system with enhanced BH 
evaporation, as it would be observed by a gravity wave detector such as the proposed Laser 
Interferometer Space Antenna (LISA). As we will discuss below, BH mass loss can lead to an 
outspiral of the binary components, while the effect of gravitational radiation produces inspiral. 
As discussed by McWilliams, for expected binary system parameters a changing orbital period 
(i.e., chirping of the gravity wave signal) could not be measured by LISA, thus one could 
not directly confirm an outspiral behavior. However, if one assumes any observed binary of 
sufficiently small orbital period has reached that period by inspiraling, then one can set a limit 
on L. For expected periods and masses, McWilliams shows that a limit of L < 5/im can be set. 

3 A Black Hole-Pulsar Binary System in the Presence 
of a Warped Extra Spatial Dimension 

We consider a binary system consisting of a BH and an NS, where the NS is a pulsar. Ob- 
servations of the pulsar could be used to measure the changing orbital period of the system 
with sufficiently high precision to either directly observe the outspiral behavior, or measure the 
competing contributions of mass loss and gravitational radiation as they affect the inspiral rate 
of the system. Alternatively, improved limits could be set on the size of an extra dimension. 
We are motivated by the PSR B1913+16 Binary Pulsar case, where observations of the one 
NS that acts as a pulsar have yielded high precision determinations of the parameters of the 
system, and have provided a dramatic test of relativistic physics. 

The BH evaporation rate in the braneworld scenario is given by 



where Mbh is the black hole mass in solar masses, and Lxq is the AdS radius in units of 10 /zm 
[7J |6] . Hereafter, M and m represent masses in units of solar masses and SI units, respectively. 
The current bound on the AdS radius is L < 44/im [T5| [T] from torsion-balance experiments. 
We take L = 10/im (i.e., L w = 1) as a nominal value in the discussions that follow. 

3.1 The Effects of Mass Loss Due to Enhanced BH Evaporation 

Consider a BH-NS binary system, with mass loss alone. We will use the results of Had- 
jidemetriou [10l ITT] , who analyzed the dynamical behavior of a binary system with isotropic 
mass loss from one or both components, in the classical case (relativistic corrections will be 
considerably smaller in magnitude than the results obtained here). Intuitively, the binary pair 
becomes less tightly bound, so the components must separate over time and the orbital period 



M BH = -2.8 x lO' 7 M BH L\ M & y 
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will increase. In a fixed, non- rotating frame of reference, let m\, and m 2 be the masses of the 
two components, and m = mi + m 2 . The rates of change of the osculating elements for the 
orbit of m 2 relative to ni\ are given by 



1 + 2ecos6' + e 2 m m 

-a = -a—, (2) 

1 — e l m m 



Tft 

e = -(e + cos^) — = 0, (3) 
m 

sin 9 m 

Co = = (4) 

e m 

where a is the semi-major axis, e is the eccentricity, u is the longitude of periastron, and 9 
is the true anomaly [TOl |TT]. The results given above are averages over one orbit, assuming 
mP <C m (as is our case); one orbit is well approximated by a Keplerian orbit where the average 
of cos 9 = — e [30] and the average of sin 9 = 0. 

Kepler's third law relates P, a, and m. Thus, to lowest order the rate of change of the 
orbital period due to mass loss (ML), averaged over one orbit, is 

P ML = 3 -P a -(l- i*?) = 2P a - (5) 
2 a \ 6m a J a 

Summarizing the results in terms of the parameters m, m, and P, we have &ml = ^ml = 
and 

G \ 1/3 / p\ 2/3 



hML = -\A^) \m) (6) 

P ML = -IP-. (7) 

m 

For mass loss, a and P increase. In the braneworld scenario the mass loss rate for the binary 
is the mass loss rate from the black hole, given by eq. ([!]). 



3.2 The Effects of Energy Loss Due to Gravitational Radiation 



The effects of energy loss by gravitational radiation from a binary system are well known 
a and P decrease — opposite results from the mass-loss scenario. For two point masses (i.e., 
a ^> radius of either object), and where averages over one orbit are of sufficient accuracy, from 
j, and [32] (in the context of the PSR B1913+16), we have 

2567r 2 G 2 m 1 m 2 (l + ie 2 + |e 4 ) 
a G R = - ' 



GR 



5 c 5 P 2 (l-e 2 ) 7/2 ' 

192tt (27rG) 5/3 m x m 2 (l + ffe 2 + §e 4 ) 
~5 8> mVsps/s T\ - e 2 ) 7/2 



(9) 
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* Gfl = 67r(4* 2 ) 1/3 V^(l-e 2 ) ("I 



Henceforth, we will concentrate on d, which is most readily observed as a rate of change of 
the period P. 



3.3 The Results for Typical Parameters: Inspiral versus Outspiral 

It is instructive to examine the above results in terms of specific values for the parameters. 
Consider a binary comprising a black hole of mass M BB - — 3 (units of M ) and neutron star of 
mass M^s — 1-4, with orbital period Piq = 1 in units of 10 hours. The semi-major axis is 

a = 2.68 x 10 9 m (M BH + M NS ) P 10 (12) 

which is ~ R Q ^> radius of either object. 

Considering the effect of mass loss only, we find 

_! / m bh + M NS y 2/3 (M BH y 2 p2/3r2 

a ML = 19 my I — 1 I — — 1 P 10 ' L 10 , (13) 

• / MBH + j^yg V 1 (M BH y 2 2 
P A/L = 0.51 msy I — I I — ^— 1 Pio^ , (14) 

which are independent of e. 

Considering only the effect of gravitational radiation, with no mass loss, we have 



-i(M BH \ (M ns \ p _ 2 , ( 
-0.38 m y I — — J I 1 Pio /(e) 



d Gfl = -0.38 m y" 1 ( -F- ) ( ^ ) Pfo7(e), (15) 
P G , = -0.0076 ms y- (*<§*) ( ( Mm±^l) ~* P ^ /(c)> (16) 



where 



1.4 / V 4.4 



/( e )=(l + Se 2 + Se 4 l(l-e 2 )- 7/2 . (17) 
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Clearly, for systems with these parameters the enhanced mass loss dominates the evolution of 
P. 

Consider a specific case with P = 7.75 hours and e = 0.6, motivated by PSR B1913+16. 
The observed rate of change of the orbital period for the binary pulsar PSR B1913+16 is 
P = (-2.4184 ± 0.0009) x 10~ 12 s/s = (-0.076 ± 0.00003) ms y" 1 [33]. For L 10 = 1, we find 
&ml = 16 m y _1 , and Pml = 0.40 ms y _1 , while the effects due to gravitational radiation are 
of opposite sign, but roughly 4 times smaller, at clgr = —3.9 m y -1 , and Pgr = —0.12 ms y _1 . 
Pml = 0.40 ms y _1 is about 13,000 times larger than the precision on the measured rate of 
change of the orbital period of PSR B1913+16. The effect of mass loss could be well-determined 
by observations of a similar precision. 
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In a BH-NS system gravitational radiation and BH mass-loss produce opposite effects. 
For small P, gravitational radiation dominates and P will decrease as the binary components 
inspiral and eventually coalesce. For large P, BH mass-loss dominates, thus P increases and 
the components outspiral. Equating the magnitudes of the rate of change of the period for 
mass loss only (eq. [T4"j) and for gravitational radiation only (eq. H~6]) . we obtain an approximate 
critical orbital period P cr iticai which separates these two disparate behaviors, 

p_ H 2 , hours L -3/ V(e) 3/S. (18) 

4 Discussion 

If a significant number of binaries have P > P cr iti ca i there are implications for discussions 
of the number of BH-NS binaries in the Galaxy, since such discussions assume all binary 
systems inspiral and eventually cease to exist. There would also be implications for gravity 
wave searches, which depend on the number of inspiraling compact object binaries. 

[22] investigated the theoretical formation rate and evolution of BH-NS binaries where the 
NS is a "recycled" pulsar (an increased spin rate from accretion, thus a greater pulsar "clock" 
stability, yielding higher precision observations, e.g., PSR B1913+16). They find that most 
newly formed BH-NS pairs of this sort have an orbital period of 2-7 hours, with a peak near 
3 hours, and an eccentricity e < 0.3, with a peak near 0.1. The BH masses range from 5-10 M , 
with an average of ~ 7 M & . 

Taking M BH = 7, M NS = 1.4, P = 3 hours, e = 0.1, and L = 10/xm, we obtain P ML = 
0.015 ms y _1 , and Pgr = —0.11 ms y _1 . The effect of gravitational radiation is ten times 
larger than the effect of mass loss, but precision measurements of the orbital period would 
show a deviation from the results due to gravitational radiation alone: Pml — 0.015 ms y _1 
is well above the obtainable observational precision of ±0.00003 ms y^ 1 (the precision stated 
here is specific to the binary pulsar PSR B1913+16 and will vary from case to case). Note that 
the critical period for inspiral versus outspiral is 6.5 hours for these masses, eccentricity, and 
AdS radius, which is within the range of orbital periods expected; it is as low as 4.2 hours for 
M BH = 5. 

[23] also determined the expected number of BH-NS binaries with a recycled pulsar in 
our Galaxy to be less than ~10. This number is quite small in part because of the expected 
coalescence rate of systems with periods of 3-4 hours and inspiral lifetimes of ~ 10 8 years. But 
the effect of mass loss, if present, is not accounted for, of course. BH-NS systems that don't 
contain a recycled pulsar are also possible, and [4] estimate there are some hundreds of such 
systems in total. Nevertheless, the search for BH-NS systems (with an observable pulsar) will 
be difficult, of course. However, searches for pulsars in binary systems are of particular interest, 
and these searches will become more sensitive with the appearance of new instruments, (e.g., 
the Square Kilometer Array). 

Given the prospects for improvement of torsion-balance measurements over the coming 



7 



decades, it is reasonable to ascertain the prospects for BH-NS observations in setting limits on 
the size of an extra spatial dimension. Taking ±0.00003 ms y _1 to be the nominal attainable 
precision for measurements of the rate of change of orbital period for any BH-NS system, then 
from eq. ([HI) , observations of a BH-NS system with those parameters could be used to set 
a 95% confidence upper limit to the size of the extra dimension of L < 0.056/im, about 800 
times better than the current 95% confidence upper limit of L < 44/im from torsion-balance 
experiments, and probably considerably better than could be set by earth-based torsion-balance 
gravity experiments in the foreseeable future. 
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